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We investigate the evolution of hydromagnetic perturbations in a small section of accretion disks. 

It is known that molecular viscosity is negligible in accretion disks. Hence, it has been argued that 
a mechanism, known as Magnetorotational Instability (MRI), is responsible for transporting matter 
in the presence of weak magnetic field. However, there are some shortcomings, which question 
effectiveness of MRI. Now the question arises, whether other hydromagnetic effects, e.g. transient 
growth (TG), can play important role to bring nonlinearity in the system, even at weak magnetic 
fields. Otherwise, whether MRI or TG, which is primarily responsible to reveal nonlinearity to make 
the flow turbulent? Our results prove explicitly that the flows with high Reynolds number {Re), 
which is the case of realistic astrophysical accretion disks, exhibit nonlinearity by TG of perturbation 
modes faster than that by modes producing MRI. For a fixed wavevector, MRI dominates over 
transient effects, only at low Re, lower than its value expected to be in astrophysical accretion disks, 
and low magnetic fields. This seriously questions (overall) suasiveness of MRI in astrophysical 
accretion disks. 


PACS: 98.62.Mw, 95.30.Qd, 47.27.T-, 47.35.Tv 


I. INTRODUCTION 

Accretion disks are found in active galactic nuclei (AGNs), around compact stellar objects in a binary system, around 
newly formed stars etc. (see, e.g., Q). However, the working principle of accretion disks still remains enigmatic to us. 
Due to its inadequacy of molecular viscosity, turbulent viscosity has been proposed to explain the transport of matter 
towards the central object. This idea is particularly attractive because of its high Reynolds number (see, e.g., 

i)- However, the Keplerian disks, which are relevant to many astrophysical applications, are remarkably Rayleigh 
stable. Therefore, linear perturbation cannot induce the onset of turbulence, and consequently cannot provide enough 
viscosity to transport matter inwards therein. 

With the application of Magnetorotational Instability (MRI; [1,0]) to Keplerian disks, Balbus & Hawley [Ij showed 
that initial seed, weak magnetic field can lead to the velocity and magnetic field perturbations growing exponentially. 
Within a few rotation times, such exponential growth could reveal the onset of turbulence. Since then, MRI has been 
a widely accepted mechanism to explain origin of instability and hence transport of matter in accretion disks. Note 
that for flows having strong magnetic helds, where the magnetic field is tightly coupled with the flow, MRI is not 
expected to work. Hence, it is very clear that the MRI is bounded in a small regime of parameter values when field 
is weak. 

It has been well established by several works that transient growth (TG) can reveal nonlinearity and transition to 
turbulence at sub-critical Re (e.g. iffil)- Such sub-critical transition to turbulence was invoked to explain colder 
purely hydrodynamic accretion flows e.g. quiescent cataclysmic variables, in proto-planetary and star-forming disks, 
the outer region of disks in active galactic nuclei. Baroclinic instability is another plausible source for vigorous 
turbulence in colder accretion disks [l^. Note that while hotter flows are expected to be ionized enough to produce 
weak magnetic fields therein and subsequent MRI, colder flows may remain to be practically neutral in charge and 
hence any instability and turbulence therein must be hydrodynamic. However, in the absence of magnetic effects, the 
Goriolis force does not allow any significant TG in accretion disks in three dimensions, independent of Re (see Q), 
while in pure two dimensions TG could be large at large Re- However, a pure two-dimensional flow is a very idealistic 
case. Nevertheless, in the presence of magnetic field, even in three dimensions, TG could be very large (Coriolis effects 
could not suppress the growth). Hence, in a real three-dimensional flow, it is very important to explore magnetic TG. 

In the present paper, we explore the relative strengths of MRI and TG in magnetized accretion flows, in order 
to explain the generic origin of nonlinearity and plausible turbulence therein. By TG we precisely mean the short- 
time scale growth due to shearing perturbation waves, producing a peak followed by a dip. By MRI we mean the 
exponential growth by static perturbation waves. While TG may reveal nonlinearity in the system, depending on Re, 
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amplitude of initial perturbation and its wavevector and background rotational profile of the flow, question is, can its 
growth rate be fast enough to compete with that of MRI? On the other hand, is there any limitation of MRI, apart 
from the fact that MRI does not work at strong magnetic fields? Note that some limitations of MRI were already 
discussed by previous authors PliMll, which then question the origin of viscosity in accretion disks. 

We show below that the three-dimensional TG dominates over the growth due to MRI modes at large i?e, bringing 
nonlinearity in the flows. This is of immense interest, as the larger Re is more plausible in accretion disks. By 
comparing modes corresponding to static (original MRI) and shearing (TG) waves, the growth estimates from static 
MRI waves have already been argued to be misleading El- Throughout their work, previous authors Ea argued 
that shearing wave structures always grow faster over short time scales than static structures, what we also plan to 
elaborate here. Nevertheless, those authors E3] did not explore the length of time over which the short-time growth 
can persist, which is very important for revealing non-linearity, which we plan to explicitly explore here. We will show 
below that in a shorter time-scale, TG reveals nonlinearity into the system. 

We explicitly calculate the magnetic field strength above which MRI not working. Moreover, for a fixed perturbation 
(which might not be either corresponding to the best MRI or best TG mode) with finite Re, with the increase of Re, 
we show that the TG tends to bring the nonlinearity in the systems before MRI could do the same, producing a large 
growth of perturbation. We notice that above a threshold Re, only TG is sufficient to make the system nonlinear 
at low magnetic held and there is no growth at high magnetic helds. Hence, in the regimes of high magnetic held 
or/and high Re, MRI is not important at all. The working regime of MRI is rather much narrower than it is generally 
thought off. As TG was argued to be plausible source of nonlinearity in cold disks and the growth due to MRI is 
subdominant compared to TG at high Re in hot disks, TG could be argued to be the source of nonlinearity and 
plausible turbulence and subsequent viscosity, in any accretion disk. 

In the next section, we discuss the perturbation equations describing hows. Subsequently, we explore total energy 
growths of perturbations due to TG and MRI for different parameter values and, furthermore, compare the respective 
parameter spaces for different initial amplitudes of perturbations in §III and §IV respectively. Finally we end with a 
discussion in SV. 


II. GOVERNING EQUATIONS DESGRIBING MAGNETIZED ROTATING SHEAR ELOWS IN 

LAGRANGIAN COORDINATES 

Within a local shearing box, in Lagrangian coordinate, the Navier-Stokes, continuity, magnetic induction equations 
and solenoidal condition (for magnetic held) can be written as 


V = —CgVp -I- -I- 2v X Q -I- -r^B • VB, 

p A-Kp 


dp 

dt 


= -pV ■ V, 


when 


^•R 

— = Vx (vxB), V-B = 0, 
dt ^ ’ 


Fiv^ 

r = v(r-), V.—.V^ 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 


where v is the velocity vector, B the magnetic held, v the kinematic coefficient of viscosity, Cg the sound speed in 
the shearing box, Q the angular velocity, p the density, r and are the position vectors in Eulerian and Lagrangian 
coordinates respectively Q. Note that the contribution of magnetic pressure has been included to the total pressure 
in the hrst term in the right hand side of equation O- For incompressible how, equation Q becomes 


V-v = 0. 


( 5 ) 


Let us dehne the tensor Qq which is the minus of the gradient of the unperturbed (background) velocity held vq = 
(0, —qflx, 0) as 


0 qn 0 

Qq = -Vvo = -(VQ) X R = I 0 0 0 

0 0 0 


9 = - 


dlnQ 

d\nR’ 


( 6 ) 
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where R = {R, 0, 0) and R is the distance of the comoving shearing box from the center of the disk and |r2| = oc i? 
(see 0 for details). Now integrating equation (|1]), we obtain 


and this gives rise to the relation 


r + ntr.q 


dr^ 

dr 


1 + ntq, 


(l + mq).V^. 


(7) 

( 8 ) 


Since the unperturbed velocity vq has spatial dependence, it has a nonvanishing time derivative in perturbed La- 
grangian coordinate. Therefore, we obtain 

dv = V — Vq = V — v.Vvo = V + riv.q. (9) 

Perturbing and linearizing equations O, ©, © and using equation ©, we obtain the perturbed Navier-Stokes, 
continuity, induction equations and solenoidal equation for magnetic field in Lagrangian coordinate as 

dv = --ctVSp + + 2(5v x Cl -\—• VdB + QSv ■ q, (10) 

p A-Kp 


Sp = —pV ■ dv, 


( 11 ) 


dB = Vx (vx(5B + dvxB) + (v-V)dB, V • 5B = 0, (12) 

where Sv, dB and Sp are the velocity, magnetic field vectors and the density of perturbation respectively. 

We now work with the incompressible approximation, i.e. Sp ^ 0 and oo, assuming c^Sp to be finite and 

decomposing the general linear perturbations into a plane wave form as 

dv, (5B oc ea:p(ik^ • r^), (13) 


when 


k = {K,ky,k^) = (1 + lltq) • k^ = (fc^ + q^tky , ky , k^), (14) 

where k and k^ are the wavevectors in the Eulerian and Lagrangian coordinates respectively. Now using solenoidal 
condition for magnetic field, incompressibility condition and plane wave solution (11311 . and if we write equations m 
and m (i.e. Navier-Stokes and magnetic induction equations) componentwise, we obtain 


Svx = —i-CgSp{k!^ 
P 


+ qfltky) — vk'^Svx + 2nSvy + -^—iSBx{Bikx + B 2 ky + B^kz) 


(15) 


Svy = -i-c^gSpky — vk^Svy — 2U,Svx + VtqSvx + - — iSBy{Bikx + B2ky + B^kz) 


(16) 


Svz = -i-clSpk^ - vk'^Svz + -^iSBz{Bikx + B2ky + Bafc^) 

P iTip 


(17) 


SBx = iSvxiBikx + B2ky + B^kz), (18) 

SBy = iSvy{Bikx + B 2 ky + B^kz) — qflSBx, (19) 

SBz = iSvz{Bikx + B2ky + Bskz). (20) 

For the convenience of solutions, we further define 

A = kxSvx + kySvy, C = kxSvy - kySvx, Ab = kxSBx + kySBy, (b = kxSBy - kySBx, 
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and for the plane wave solutions given by equation 0, equations (Ull)-® can be recast into 


where 


1 ^ ^ 

/ 

c 


A_b 


\ Cb / 

V 


Mil Mi 2 Mi 3 iVfl 4 
M 21 M 22 M 23 M 24 

M 31 M 32 -M 33 M 34 



Mil = 
Mi 3 = 


—vk^{k'^ + ky) + 2qVlkxkyk‘l 

P(/c2+fc2) 

i{Bikx + B2ky + B^k^) 


Ml 9 = 


2»fc2|fe2 + (1 _ g)fc2| 

fc2(fc2 +fc2) 


47rp 


Mi4 — 0, 


M 2 l=n{q- 2 ), M 22 = -vk\ M23=0, + 

47rp 

-^31 = i{Bikx + B 2 ky + B^kz), M 32 = 0 , M 33 = 0 , M 34 = 0 , 


qQ(k'^ — k"^) 

Mil = 0, M42 = i{Bikx + B 2 ky + B^kz), M43 = — ^ — 

+ ^ 7 ; 


Mii — 


2qQkxky 

kl + kl ■ 


( 21 ) 


The assumption of incompressibility is justified as follows. If the wavelength of the velocity perturbations is much 
shorter than the sound horizon for the time of interest (which is in the present context the infall time of matter), then 
the density perturbations (which is basically the sound waves) reach equilibrium early on, which renders effectively 
a uniform density during the timescale of interest. For an astrophysical accretion disk around a black hole, which 
is either geometrically thin or can be approximated as a vertically averaged flow, the half-thickness of the disk is 
comparable to the sound horizon corresponding to one disk rotation time. Therefore, as described in previous work 
(e.g. 0), for processes taking longer than one rotation time, wavelengths shorter than the disk thickness can be 
approximately treated as incompressible. 

Solving the set of differential equations (I2ip . we can calculate i5v, dB and the energy £ of the perturbation given 

by 


£ oc 



/47rp(A2 + C2)+A|+C| 4^A2p+A|\ 

{kx + kyqmf + kl ^ kl 


( 22 ) 


in terms of new variables. In order to solve the set of equations m, we have to supply (5B and Sv at t = 0, i.e. initial 
perturbation amplitude (IPA). The structure (and evolution) of perturbations are similar/same as that found earlier 
E01I3- A sample is shown in Fig. [TJ demonstrating how an initial leading wave, with a highly stretched structure, 
evolves to a spherical wave at the maximum of TG and furthermore evolves to a trailing wave, during the declining 
phase of TG. During this evolution of perturbation, observing the associated total energy growth of perturbation, 
we now plan to understand whether the perturbation will sustain or not to give rise to nonlinearity and plausible 
turbulence and essentially viscosity to help infall of matter in an accretion disk. By a detailed investigation, we can 
also understand the relative weight between TG and growth due to MRI (if at all working) in the time of interest. 
Moreover, we plan to pinpoint the limit of magnetic field strength, above which the MRI is suppressed (indeed MRI 
works only for weak magnetic fields). 


III. TOTAL ENERGY GROWTH OF PERTURBATIONS FOR DIFFERENT PARAMETER VALUES 

The best possible mode for MRI giving rise to the nonlinearity into the system corresponds to the condition 
kzVAzI^ = I (when v\^ = Bl/Airp) 0. The growth rate for this fastest exponentially growing mode is 3D/4 = 3/4q 
(since in dimensionless unit = I/g) 0il3. Is there any mode for which TG brings in the nonlinearity into 
the flow (the best possible mode for TG) at a timescale shorter than the rotational time at which the best possible 
MRI mode brings in the nonlinearity? Note that an approximate emergence of nonlinearity is defined through the 
measurement 


IH , W 

ivi |B|; ■ 


Linearity 


(23) 
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FIG. 1: Development of the perturbed velocity 5vx{x,y) as a function of time, when Re — 10®, ky = I, z = 0 and tmax 
denotes the time at which growth attains its maximum, at (a) t = 0, (b) (c) tmax and (d) Ztmaxl2. The gradual 

conversion of contour colors from white to black corresponds to the gradual conversion from positive to negative values of 
5vx{x,y) respectively. 



FIG. 2: (Color online) Nonlinearity via best possible TG and MRI. Thick black line corresponds to the TG for IPA= 10“®, 
Re — — —RI^^, ky = 1, kx = dotdashed black line corresponds to the TG for IPA= 10~®, Re = 10^®, 

k^ = —ky = 1, kz = 90fc,^; red longdashed and dotted lines correspond to the best possible MRI starting from IPA 
= 10“® and 10“® respectively. Dashed horizontal line indicates linearity unity. 
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FIG. 3: Same as Fig. but the black thick and dotdashed lines correspond to the TG for = 1, ky = 1, kz = 100, Re = 10^^. 
and k^ = 1, ky = 1, kz = 3000, Re = 10^^ respectively. 


When Linearity=l, the system will start becoming nonlinear which will plausibly lead to turbulence. For a Keplerian 
disk {q = 3/2), the best MRI mode brings in the nonlinearity at the timescales ~ 14 and 23 respectively for IPAs 
= 10“^ and 10“® (when log(l/IPA) = it/Aq). Figure[2]shows that indeed there are modes which reveal nonlinearity via 
TG at around 3 and 13 rotational times for IPAs 10“^ and 10“® respectively (where |(5v(0)|/|v|, |(5B(0)|/|B| = IPA). 
Note from Q that the maximum TG in two dimensions scales as TGmax and the corresponding time as 

in pure hydrodynamical disks. This further corresponds to k^: ^ [l^, which reveals TGmax corresponds to the 

minimum of .mzn in. In the same spirit, k^s in Fig. [2] are chosen to be when initial perturbations 

are highly stretched and nonspherical. Note that although such a stretched initial wave vector of perturbation is a 
special choice which is important for the present purpose, nothing prevents them to be arising in nature. Since every 
perturbation mode is equally probable when a system is perturbed (which is indeed the idea behind the choice of the 
best MRI mode), we explore the mode which is growing faster and leading the system to nonlinearity. In Fig. |3l we 
relax this Re dependence of k^, but still obtain the nonlinearity arising at ^ 3.5 and ^ 6 rotational times for IPAs 
10“^ and 10“^ respectively. Hence, the full-scale general hydromagnetic effects giving rise to TG reveal nonlinearity 
into the system faster than that the MRI does, when MRI itself is uncertain. Once the best TG reveals nonlinearity 
before the best MRI would do, the importance of MRI is sluggish in the linear theory. Note that our current emphasis, 
in particular, is the arisal of nonlinearity via either TG or MRI. However, nonlinearity does not guarantee for the 
transition to turbulence, in the physical time scale of accretion. One could argue that MRI modes grow forever and, 
hence, the system would have been turbulent at some point, even if the best MRI modes are not considered when 
TG would eventually decay. But the important fact to notice here is that as soon as the system becomes nonlinear, 
MRI (and also TG, if that wins over MRI) is no longer applicable, as the underlying solution itself is based on linear 
theory. On the other hand, the effects due to the best possible MRI mode should be compared with that of the best 
possible TG. Such a comparison shows that TG is more powerful and is actually responsible for bringing nonlinearity 
into the systems. 

Let us move on to the detailed behaviors of TG. In Fig. 01 we show energy growths {£{t)/£{0)) for four sets of Re 
and (plotted in dimensionless units, based on the dimensions of various quantities of shearing box). We see that 

for a fixed Re-, energy growth of perturbation decays over time if the background magnetic field is large (thick and 
longdashed lines compared to the dotted and dotdashed lines respectively). Figure [5] shows the linearity of respective 
cases based on equation (l23l) . The most important point to be noted from Fig. [5] is that the case of high Re and low 
B (dotted line) exhibits nonlinearity via TG itself, for IPA ss 10“^. Note the clear appearance of TG peak in the 
linearity as well as growth curves at time t ^ 10^. Later, at the trailing phase of TG, growth further starts increasing 
due to MRI. However, by this time, the system would already become nonlinear, and hence computations of energy 
growth and that due to MRI based on the linear theory lose their meaning. However, for a lower Re for the same B, 
growth due to MRI overpowers TG and nonlinearity arises via MRI induced growth (dotdashed line). 

In order to understand a global picture and the relative powers of TG and MRI, we perform several numerical 
experiments and in Fig. [Slwe divide the entire parameter space into three regions: MRI active, TG active and stable 
(or linear) zones, for a given perturbation wavevector. Note that for a given B, the difference of log(i?e) between two 
successive computations is chosen to be unity and hence the curve dividing the linearly stable (no energy growth of 
perturbations) and unstable zones does not appear very smooth. The region left to the solid vertical line exhibits 
nonlinearity via MRI, while that of right side corresponds to nonlinearity via TG before MRI could kick in. 

Hence, a very important message here is that energy growth rate due to MRI is faster than TG only at lower values 








7 



FIG. 4: Total energy growth for different sets of Re and B = (0,0,53): Thick, longdashed, dotted and dotdashed lines 
correspond to respectively Re = 10^^ and /p = 10~®, Re = 10^ and /p = 10, Re = 10^^ and B^/p — 10”^'^ and 
Re = 10^ and B^/p = 10“^°. = —RI^^, ky = kz = 1. Inset confirms that the oscillatory zone of thick line is continuous and 

smooth. 


c 

3 



FIG. 5: Linearity of cases in Fig. [4] Dashed horizontal line indicates linearity unity. Inset confirms that the oscillatory zone 
of thick line is continuous and smooth. 


of Re and it is further suppressed above a certain higher B (when indeed MRI is a weak field effect). At larger Re, 
which actually corresponds to astrophysical accretion disks, growth rate due to TG overpowers that due to MRI. 
Although Fig. |6] represents cases corresponding to a vertical background magnetic field, we obtain similar trends of 
results at other background magnetic field geometries and wavevectors, and hence they are not shown here. 


IV. COMPARISON OF PARAMETER SPACE EOR DIEEERENT INITIAL AMPLITUDE OE 

PERTURBATIONS 

Figure [7] compares the parameter spaces, as described in Fig. |6l for two different IPAs: 10“^ and 10“^. As IPA 
decreases, the value of Re dividing the MRI active and TG active zones, namely Red, increases, which apparently 
implies that the MRI active region increases. It also appears that Red ~ IPA“^. Let us now recall the time scale 
leading to the system to be nonlinear by the respective growths due to MRI and TG and estimate if those tally 
with observation and initial choice of our model. For this purpose, first we fix Re at 10^^ and take the IPA to be 
10“^, along with sufficiently low B so that the flow is assured to be nonlinear and unstable, in the parameter space 
described in Figs. [5] and [T] In this case, nonlinearity arises in the TG active zone at about 750 rotation time and 
corresponding TG is shown in Fig. [5] (dotted line). If the shearing box is at lOORg away from a IOMq black hole, 
where Rg and Mq are Schwarzschild radius and Solar mass respectively, then this dimensionless time scale recasts 
into 750L/{qil,L) = 750y/R^/GMq'^ ^ 750 seconds for q = 1.5 (Keplerian disk), when L is the radial width of the 
shearing box, G the Newton’s gravitation constant and M the mass of the black hole. Now if we decrease IPA to 
10“^ keeping Re fixed, TG cannot bring nonlinearity any more (however by increasing Re, again TG could bring 
nonlinearity), as shown in Fig. [71 instead, the nonlinearity arises then via MRI modes. However, the time scale 
























FIG. 6: Parameter space describing stable and unstable zones, based on the MRI and TG inactive and active regions, when 
B = (0,0, Bs) and IPA is 10“®. The dashed vertical line at Re = 10® is the threshold Re above which MRI does not work. 
ky = kz = 1 . 



FIG. 7: Same as Fig. [6] but comparing results with different IPA, when solid and longdashed lines are for IPA = 10 ® and 
10“® respectively. The dashed and dotted vertical lines at Re — 10® and 10^® correspond to threshold Re (as in Fig. O when 
IPA = 10“® and 10“® respectively. 


for the emergence of nonlinearity, in this MRI active process, as shown in Fig. [5l is approximately 35000 rotation 
time which is ^ 35000 seconds (following the same procedure, as used above for the TG active case, to convert the 
dimensionless to dimension-full times). Now we can calculate the radial velocity (vr) of the Keplerian accretion disk 
at the location of shearing box for a given accretion rate M, say 0.1 Eddington rate [l^, which is supported by 
observation, given by 


Vr 


2 X 


/ M / M ( R 

1^3 X 10-8Mo/year ) [m^J [jlfg) 



(24) 


when a is the Shakura-Sunyaev viscosity parameter (whose origin is actually aimed here to determine). It is reasonable 
to assume that the time required to make the flow nonlinear, and hence turbulent which subsequently reveals viscosity, 
is of the same order as the time required by a fluid parcel to cross the length of the shearing box radially (t^) as a result 
of turbulent viscosity. Hence, the product of Vr and should be of the order of width of the shearing box L. For the 
above mentioned case of IPA= 10“^, when nonlinearity is due to TG, we obtain L = 0.1i?g from equation (1^^ . which 
is highly reasonable for our choice of shearing box approximation (L « R). However, for the case of IPA= 10“®, 
when nonlinearity is due to MRI, we obtain L = lORg which marginally satishes (or even violates) the initial choice 
of a narrow shearing box at lOORg. Therefore, although smaller IPAs increase MRI active zones, the observed infall 
cannot be explained by them. This problem with MRI would appear to be more severe at progressively lower IPAs 
and TG would be more important for revealing nonlinearity at progressively higher Rg, which are forbidden for MRI. 
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V. DISCUSSION AND CONCLUSIONS 

Let us estimate the maximum |B| in Gauss supporting nonlinearity, as shown by the solid curve(s) in Figs. |6]and 
[71 We again set the shearing box at 100i?g away from a IOMq black hole. Then we obtain the values of density 
(PiooRg) at that location to be 10“^ gm/cc [^. The background Keplerian velocity at that position, for the 
size of the shearing box, 0.1i?g, which is consistent with that obtained for the TG active zone, can be obtained as 
qflL = GM/R^L ~ 10® cm/sec. We now consider = 10^^ and, hence, from Fig. |6]the corresponding maximum 
(dimensionless) magnetic field supporting nonlinearity is given by B'^jp = 10“®. Therefore, corresponding actual 
value of magnetic field is piQOR^{q^lLY ^ 30 Gauss. This means, the flow with = 10^^ and |B| >30 Gauss, 

the energy growth of perturbation will decay over time, but for |B| ^ 30 Gauss, TG will be sufficient enough to bring 
nonlinearity in the system, however, still not requiring any growth due to MRI. ^From Fig. |6l it is clear that MRI 
is only important whenever Rf. < 10®, whereas for Re ^ 10®, which is the favorable zone of Re for accretion disks, 
magnetic TG is important than MRI. 

In short, we have calculated the magnetic field strengths for different ReS above which the system will be stable 
under linear perturbation and an upper bound of Re above which either the system is stable under linear perturbation 
(for high magnetic field strength) or reaches nonlinear regime (for low magnetic field) through magnetic TG. In one 
line, MRI is not the sole mechanism to make accretion disk unstable, there is a big kingdom where TG rules, and 
explanation of accretion solely via MRI is misleading. 

The authors would like to thank Prateek Sharma for illuminating discussion. B.M. acknowledges partial support 
through the research grant provided by Indian Space Research Organization of Ref. No. ISRO/RES/2/367/10-11. 
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